Competitive Clustering in a Bi-disperse Granular Gas 
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A bi-disperse granular gas in a compartmentalized system is experimentally found to cluster 
competitively: Depending on the shaking strength, the clustering can be directed either towards 
the compartment initially containing mainly small particles, or to the one containing mainly large 
particles. The experimental observations are quantitatively explained within a flux model. 

PACS numbers; 45.70.-n 



Clustering is one of the important features of granular 
gases, arising from the inelastic collisions between the 
particles. Energy is dissipated in each collision such that 
a dense region will dissipate more energy, becoming even 
denser, resulting in the formation of a cluster of slow 
particles Q . It is a prime example of structure formation 
in a system far from equilibrium. 

Thus far, most attention has been given to clustering in 
mono-disperse systems where all particles are identical. 
In that case, the clustering simply occurs in a region 
that originally is a bit denser than the others. In the 
present paper we concentrate on a bi-disperse granular 
gas, consisting of large and small particles. Here the 
clustering turns out to be competitive: it can occur either 
in a region which originally is populated mainly by small 
particles, or in a region with mainly large particles. 

The experimental setup (see Fig.|l|) consists of a cylin- 
drical perspex tube with inner diameter 11.2 cm and 
height 42.2 cm, divided into two equal compartments 
by a wall of height 6.0 cm. This compartmentalizaton 
makes it possible to get a clear-cut picture of the clus- 
tering process. The tube is mounted on a shaker with 
adjustable amplitude a and frequency /. The inverse 
shaking strength D oc \/{af)'^ (defined more precisely 
in Eq.(|8|) below) is one of the crucial parameters for the 
clustering behavior. 

To minimize the effect of statistical fluctuations in the 
experiments, we take a sufficiently large number of beads 
of each species: Pi = 300 large and P2 = 600 small ones. 
The small and large beads are both made of steel (so they 
have the same coefficient of normal restitution, e w 0.85), 
with radius ratio = rxjri = 2. 

Clustering behavior — The clustering behavior is 
strongly influenced by the initial distribution of the two 
species over the compartments. As an example we take 
an initial distribution in which compartment A has the 
majority of the large particles, and compartment B most 
of the small ones. This situation is depicted in Figj^, 
left, with {|Pi,iP2} = {180, 200} in compartment A 
and {|Pi,|P2} = {120, 400} in compartment B. The 
outcome of the experiment (see Fig.|l|, right) depends on 
the shaking strength. 

The full bifurcation diagram is given in Fig.||: For suf- 
ficiently strong shaking (small D) a uniform distribution 
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FIG. 1: Images from two experiments with a mixture of 
Pi =300 large steel beads (ri = 2.50 mm, mi — 0.500 g) 
and P2=600 smaller ones (r2 = 1.25 mm, m2 = 0.063 g). The 
initial distribution of the beads is the same in both experi- 
ments: {|Pi, IP2} in compartment A and {f Pi, IP2} in B. 
For moderate shaking strength (/ = 60.0 Hz, a = 2.0 mm, top 
row) a cluster is formed in the compartment initially domi- 
nated by large particles (A). For milder shaking (/ = 37.5 Hz, 
a — 2.0 mm, bottom row) the clustering takes place in the 
compartment initially containing most of the small particles 
(B). 



is found. This is denoted as regime 0. For moderate 
shaking strength (moderate D, regime I) the particles 
cluster in compartment A, the one initially dominated 
by large particles. The reason is as follows: Many of the 
small beads quickly cluster in compartment A, where the 
dissipation is highest due to the greater number of large 
beads, which (with their larger mass and surface area) 
act as "coolers" . The remaining beads in compartment 
B jump higher than before, since there are fewer colli- 
sions, and also large ones now start to make it over the 
wall into compartment A. After about twenty seconds to 
a minute (depending on the value of D) the final state 
is reached: a dynamical equilibrium with practically all 
large particles and most of the small ones in compartment 
A, and only a few rapid small particles in compartment 
B (see FigJ^, top row). 
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FIG. 2; Bifurcation diagram showing the three different clus- 
tering regimes O, I, and II. The particle numbers Ni (i — 1, 2) 
are given relative to the symmetric solution: Ni — \Pi- The 
squares and diamonds are experimental data, corresponding 
to the numbers of large and small particles respectively, in 
compartment A (solid) and B (open). The curves represent 
the equilibrium situation according to the theoretical flux 
model. The initial condition is always the same as in Fig|l[ 



For very mild shaking (high D, regime II) the clus- 
tering process is much slower and, more importantly, it 
goes in the opposite direction. The particles now clus- 
ter in compartment B, i.e., in the compartment initially 
dominated by small particles (see Fig.|l|, bottom row). 
The series of events is as follows: First, the larger par- 
ticles hardly jump at all. They stay close to the bot- 
tom, transferring energy from the vibrating bottom to 
the smaller particles, which thereby gain comparatively 
large velocities (like tennis balls jumping higher on top 
of a basketball than on the plain floor). So it is easier 
for the small beads to leave compartment A (with more 
large beads). The remaining particles become more mo- 
bile, and after a couple of minutes the first large beads 
start to make it across the wall into B, where they are 
immediately swallowed by the developing cluster. With 
every particle that leaves compartment A, the process 
progressively speeds up. The total clustering time is in 
the order of five to twenty minutes, steeply increasing 
with the value of D. 

Flux model — The above behavior can be explained 
in terms of a dynamical flux model. It is a generaliza- 
tion of the model derived by Eggers for clustering in the 
mono-disperse case H, ^ . First we consider the com- 
partments separately, and from this derive the flux of 
particles leaving each compartment. 

The gas within a compartment is assumed to be in 
thermal equilibrium with respect to the granular tem- 
perature T = (wf ) (the mean kinetic energy of the 
particles) , which is taken to be the same for both species 
i — 1,2. This is an approximation: Two different granu- 
lar temperatures may in fact coexist in binary mixtures 
B (and at higher densities the species are known to seg- 



regate p), but taking this into account would add an 
unknown and possibly not constant parameter to our 
model. As in ref. and in fair agreement with molecu- 
lar dynamics simulations , we assume that the granu- 
lar temperature T is approximately constant throughout 
the compartment. 

The same simulations |7| show that both species obey 
a barometric height distribution (as expected for dilute 
granular gases, where one can use the standard equation 
of state and momentum balance [^). That is, their num- 
ber densities show an exponential decay with the height 



ni{z) = ?ii(0) exp(-mi5z/T). 



(1) 



The density at ground level follows from Vl Jg°° ni{z)dz — 
Ni as rti(O) — niigNi/ilT. Here Ni is the number of 
particles (of species i) in the compartment under con- 
sideration and fl is its ground area. Eq.Q says that the 
larger particles, with larger rrii, become dilute faster with 
growing z than the small particles. 

The temperature T is determined by the balance be- 
tween energy input Jo (through collisions with the bot- 
tom) and dissipation Q (through the collisions between 
particles). To determine the energy input, we assume 
for simplicity a small-amplitude sawtooth motion of the 
bottom, such that colliding particles always find it mov- 
ing upward with velocity Vb = af. This means that 
a particle coming down with vertical velocity compo- 
nent Vzi is bounced back with Vzi + 2vb, and the energy 
gain per collision is 2miVb{vzi + fb). Multiplying this 
by the number of collisions (oc ni(fi)vzi per species), and 
assuming an isotropic Maxwellian velocity distribution 
{(vli) = I (vf)), this gives the rate of energy input: 
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As typically Vb ^ Vi, the first term is much larger 
than the second, which we therefore neglect. Then 
Jo = |l7(ni(0) -f n^mvbT = ^gVb{miNi + m^N^). 

This must be balanced by the rate of energy loss 
Q through the inelastic collisions between the parti- 
cles. It is equal to the product of the number of colli- 
sions per unit volume and the energy loss per collision 
(cx (1 — e'^)[raim.j I [mi + mj)][vi — fj]^), averaged over 
the velocity distribution, and integrated over the whole 
compartment. It consists of three terms, representing 
collisions between two particles of type 1, type 2, and 
type 1 and 2, respectively [§1: 
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Equating Jq and Q yields the granular temperature T 
of the compartment: 



T = 



1447r(l - e2)2 ' 
where the effective mass ^ is given by: 

n{miNi + m2N2) 



(5) 



(6) 



It is through this quantity that the populations of the two 
species within one compartment influence each other. 

The central quantity of the model is the flux function 
Fj, i.e., the number of particles (of species i) that leave 
the compartment per unit time. It is the product of the 
density 'ni[z) and the horizontal velocity (= 2T / inii) , 
integrated over the space above the wall (width h) from 
z = ft. to some cut-off height ft + iJ : 
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In the last step we have linearized exp{—migH/T), im- 
plying that H <C (vf) /g. The prefactor determining the 
absolute rate of the flux is given by ^ = 12^27r/3(l - 
e^)gbH/flaf. The dimensionless parameter D, which 
governs the clustering behavior, has the form 



D = 1447r- 



gh 



(8) 



The evolution of the number of particles NiA in com- 
partment A (j = 1, 2) is given by the net balance between 
the (outgoing) flux from A to B and the (incoming) flux 
from B to A: 



dN, 



dt 



- -F,{NiA, N2a) + F,{Nib,N2b) 
=. -F,{NiA, N2a) + F,{Pi ~ NiA, P2 " N2a), 



(9) 

where we have used particle conservation, NiA + A^is = 
Pi . The evolution of the (complementary) particle num- 
bers in compartment B is governed by the same equation 
with A and B interchanged. 



Flux model vs. experiment — In Figj^ the results from 
the flux model are shown together with the experimental 
data, using the same initial distribution as in Fig.|l]. The 
agreement between theory and experiment is very good, 
and the same three regimes are found: For vigorous shak- 
ing (regime 0, D < 10) the system settles into the unclus- 
tered homogeneous state. For moderate shaking (regime 
I, 10 < D < 150) we find clustering in compartment A, 
and for mild shaking {D > 150) the clustering takes place 
in compartment B. The transition from regime I to II is 
quite abrupt. 

The experimental observation that the small particles 
always (in both regimes I and II) are the first to clus- 
ter, followed later by the large ones, is also reflected in 
the flux model. This can be seen in Figjl, where the 
evolution of the various particle numbers (obtained from 
Eq.^) is shown. The same plots also illustrate the dif- 
ferent timescales for type-I and type-II clustering: The 
time it takes the system to reach its equilibrium situa- 
tion grows rapidly with growing D (decreasing shaking 
strength), just as in experiment. 

In the lower part of Fig.^ the corresponding flow di- 
agrams are depicted, showing how the particle numbers 
Nib and N2B in compartment B evolve, for any initial 
condition. The cross denotes the initial condition used 
in the experiments {{Nib, A''2b}={120, 400}), and the ar- 
rows show the evolution of the system. For very strong 
shaking (left plot) only one stable fixed point exists: the 
uniform distribution {150,300}. 

In the middle plot, for moderate shaking strength, the 
homogeneous state has become unstable and has given 
way to two new stable fixed points. These correspond 
to the compartment B being either nearly empty (fixed 
point in the lower left corner) or well-filled (upper right). 
Their basins of attraction are indicated by the shading. 
Our starting point happens to lie in the basin of the lower 
left point, and so we end up with compartment B being 
nearly empty. The arrows show that the small particles 
cluster first, and only when these have attained their final 
distribution do the large ones follow. Indeed, the fiux of 
large particles is only appreciable when the number of 
small particles in the compartment has become very low. 

Finally, in the rightmost plot (for mild shaking) we see 
that the basin boundary has shifted, which explains the 
switch from regime I to regime II: Our initial condition 
now lies within the basin of attraction of the fixed point 
in the upper right corner. So we end up with a full com- 
partment B. The same plot shows that the fixed points 
move further into their corners as D grows, i.e., the clus- 
tering becomes more pronounced for decreasing shaking 
strength, just as in the mono-disperse case H, 

Further parameters — Until now we have taken the 
ratio of the radii ip — ri/r2 — 2, but how does the 
phenomenon depend on ip in general? This ratio has 
a marked effect on the critical Z?- values where the transi- 
tions between the regimes 0, 1, and II take place. In Fig.0 
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FIG. 3: Evolution of the bi-diperse system, calculated from the flux model, for (a) D = 1, (b) D = 100, and (c) D = 200. 
The top rows show the number of large and small particles as a function of time (note the different time scales): the solid 
curves represent compartment A, the dashed ones compartment B. The flow diagrams (bottom row) show how the contents of 
B evolve: Nib (large) and N2B (small). The cross indicates the initial condition {A'^ib, Af2s}={120, 400}. The transition from 
(a) to (b) comes about through a bifurcation, whereas the transition from (b) to (c) is seen to be the result of a shifting basin 
boundary. 
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FIG. 4: Phase diagram, showing the three clustering regimes 
as a function of the inverse shaking strength D and the size 
ratio 1/) = ri/r2. The shaded areas are calculated from the 
flux model, and the open circles (no clustering), stars (type 
I clustering), and triangles (type II clustering) correspond to 
experiments. The experimental results on the vertical dashed 
line ip = 2 also feature in Fig.^. The initial distribution of 
particles is always taken to be the same as in Figs. pJ-H. 



we show the position of these regimes as a function of t/i, 
from the flux model and from experiments. The vertical 
dashed line corresponds to i/' = 2 studied in Figs. 1-3. 
It is seen that for i/; < 1.5 the transition to regime II is 
immediate: here the larger beads are not sufficiently big 
to compensate for the fact that they are a minority. It is 
the larger number of beads that decide where the cluster 
goes, just as for the mono-disperse case {tp = 1). On the 
other hand, for high values of ■0, the dominant size of the 
large beads always makes them the decisive factor (only 
regime I survives). It is precisely the intermediate region 
\.h < Ip < 2.3 in which the competition takes place. For 



■0 « 1.6 we witness the particularly interesting sequence 
O-II-I-II, both in the model and in experiment. 

Another parameter of interest is the ratio between the 
total particle numbers of each species, a = P1/P2 (up 
to now we have always taken a = 1/2). Also this ratio 
obviously has a big influence on the clustering behavior: 
a larger value of a means that the large beads become a 
more important minority (or even a majority for a > V) 
and hence type-I clustering will gain ground. 

Conclusion — The clustering behavior of a bi-disperse 
granular gas is much richer than the mono-disperse case, 
and shows a surprising new feature: The clustering can 
be directed either towards the compartment initially con- 
taining mainly large particles (type-I clustering) or to 
the one containing mainly small particles (type-II), sim- 
ply by adjusting the shaking strength. All the experi- 
mental observations can be explained quantitatively by a 
bi-disperse extension of the Eggers flux model. 
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